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This paper is a second in a series devoted to the study of a two- 
oscillator system in linear relative motion (the first one published as 
a letter in Europhys. Lett. 91, 60003 (2010)). The main idea behind 
considering this kind of system is to use it as a simple model for 
Casimir friction. In the present paper we extend our previous theory so 
as to obtain the change in the oscillator energy to second order in the 
perturbation, even though we employ first order perturbation theory 
only. The results agree with, and confirm, our earlier results obtained 
via different routes. The friction force is finite at finite temperatures, 
whereas in the case of two oscillators moving with constant relative 
velocity the force becomes zero at zero temperature, due to slowly 
varying coupling. 

PACS numbers: 05.40.-a, 05.20.-y, 34.20. Gj 

1 Introduction 

Consider two dielectric or metallic slabs with parallel surfaces, closely sepa- 
rated. If the slabs are set in tangential motion with respect to each other, 
with constant velocity, they become exposed to a friction force called Casimir 
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friction, an effect that has received considerable attention in the recent past. 
The physical reason for the effect is that photons transferred between the 
slabs are subject to Doppler shifts. Such frequency shifts are physically ex- 
pected to lead to energy dissipation, and hence a friction force. Some papers, 
limited to the period from 2007 onwards, are listed in refs. [TJ |2j |3l HJ El 
El El HP] . As one might expect, the problem becomes singular if the relative 
velocity of the slabs is simply taken to be constant, from t = — oo to t = +oo. 
It is physically advantageous, therefore, to imagine that the interaction is ef- 
fectively coupled in for a finite period of time. We shall make use of this 
technique in the following. It is in principle similar to the technique often 
employed in scattering field theory. 

Now it might appear most natural to attack the Casimir friction prob- 
lem by making use of standard macroscopic electrodynamics for media in 
constant rectilinear motion. It means use of Maxwell's equations in moving 
matter endowed with a refractive index n. And as expected, it turns out 
that most of the mentioned papers are following this kind of approach. A 
complicating factor in the present case is, however, that there is no obvious 
rest system of the matter to refer to; none of the slabs plays a privileged role. 
This contrasts the usual case in phenomenological electrodynamics where the 
covariant formulation is simply constructed such that the theory reduces to 
conventional electrodynamics in the uniquely defined rest inertial system (cf., 
for instance, Refs. [Ill [EH CE3] ) • 

Faced with this circumstance or difficulty, it lies at hand to search for al- 
ternative approaches to the friction problem. On natural possibility is then 
to make use of statistical mechanical methods, for harmonic oscillators in 
uniform relative motion. These methods were used by us in a recent work 
[H] , and were used also in previous works [El [TB] . The microscopic method 
has some advantages in comparison with macroscopic electrodynamics, most 
notably that the formalism becomes more simple and transparent. Yet, the 
statistical mechanical approach has the property that it is capable of show- 
ing the main features of the problem. The microscopic approach has also 
been followed recently by Barton (preprint, personal communication). In 
the present paper we follow basically our earlier avenue of research. 

To begin with, let us briefly summarize the basic characteristics of the 
model we consider (cf. also Ref. [H]): 

The two-oscillator system is quantum-mechanical in nature; the reference 
state corresponding to a Hamiltonian Hq with no coupling. Thermal equi- 
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librium is initially assumed. The equilibrium state becomes disturbed by a 
time-dependent term written as —Aq(t), where A is a time- independent op- 
erator and q(t) a classical function of time. Its explicit form depends on the 
details of the system. The Hamiltonian now becomes H = Ho — Aq(t). As 
a specific case we considered in Ref. [H] the situation with two oscillators in 
relative motion with 

- Aq{t) = ip(r) Xl x 2 , (1) 

ijj(r) being the coupling strength, r the separation between the oscillators, and 
x%,x 2 the internal vibrational coordinates. If v denotes the non-relativistic 
constant relative velocity, the interaction varies as 

- Aq(t) = [V>(r ) + W(r„) ■ vt + ..]x x x 2 , (2) 

when expanded around the initial position r = ro at t = 0. The force B 
between the oscillators is 

B = -(W)xi2; 2 . (3) 

However, as mentioned above the situation becomes singular, so to obtain 
the dissipation the perturbing interaction should last for only a finite period 
of time. Thus vt in Eq. (j2j) was replaced by \te~ vt for t > (and for 
t < 0) and the limit 77 — > was considered. For that situation we found no 
dissipation at T = 0. Then a more general interaction —Aq(t) was considered 
in Ref. [14], and an expression for the dissipation was found where friction 
could be present also at T = 0. To arrive at this result we also associated the 
general situation with motion. In view of the attention paid to the friction 
problem it is of interest to check the results of Ref. [14] in an independent 
way, as we do in the present work. 

After this brief survey let us now point out the chief results obtained in 
the present paper: 

• We calculate a general expression for the change in energy, called AE, 
by means of quantum mechanical perturbation theory. This energy 
change occurs to second order in the perturbation. Nevertheless, time- 
dependent perturbation theory to the first order is sufficient to find the 
second order effect. This is because the phases of the perturbed change 
in amplitude, and the initial amplitude of each eigenstate, are uncorre- 
cted at thermal equilibrium. Change in the amplitudes of eigenstates 
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will accordingly be the square of perturbed amplitudes; i.e. there are 
no cross-terms. We moreover find that AE is always positive or zero, 
corresponding to a friction force. Doubts occasionally raised in the lit- 
erature about the very existence of the Casimir friction effect [5] are 
thus from this standpoint laid at rest. 

• Making use of the mentioned expression for AE we compare the present 
formalism with that of Ref. [H], the latter result obtained in a quite 
different way and on a quite different form. In Ref. [H], the linear 
response via the Kubo formalism was used [T5J [TBI ELT] to calculate the 
force which in turn could be divided into a reversible and an irreversible 
part. It is the latter part that is associated with dissipation. A sat- 
isfactory feature is that the present derivation, although being quite 
different from that of Ref. [H], leads to the same physical result. 



2 Time-dependent perturbation theory 

To fix the notation, we consider the perturbation theory for a system at 
thermal equilibrium. The wave function can be written as 

j> = ^ a n ipn, (4) 

n 

where ip n = i/) n (x) are the eigenstates. For simplicity we here let x represent 
all the coordinates of the system. If ip is normalized, J ip*ipdx = 1, then 
|a n | 2 is the probability for the system to be in eigenstate n. At thermal 
equilibrium this probability is given by the Boltzmann factor 

Pn = \ a n\ = e ^ "i (5) 

where E n is the energy eigenvalue of the state and Z is the partition function 

Z = J2^ En - (6) 

n 

Let now the Hamiltonian be perturbed by the time-dependent interaction 
V(t) = —Aq(t), as mentioned above. Due to the perturbation the coefficients 
a n will change. If the system starts in a state m there are transitions to other 
states given by a change in a n 

Aa n = b nm . (7) 
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The b nm is given by the standard expression 
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h nm = -j. / V nm (r)e luJ - T dr, (8) 



OO 



where 

V nm (r) = J i)* n V{r)i) m dx = -A nm q(r), 

A nm = (n\A\m) = J ^ n Ai) m dx. (9) 

Here u nm = u n — co m , with u n = E n /h. 

As mentioned above, to avoid singularities caused by too idealized condi- 
tions we will assume that the perturbation vanishes after some time. Then 
we will obtain the total change in Aa n with 

in 

/oo 
q(t)e- iujt dt, (10) 
-oo 

where the hat denotes Fourier transform. 

From a general perspective, the system may start in a combination of 
eigenstates with transitions from several states. (It might be natural in this 
context to think about the Casimir-Polder setup with molecules traveling 
close to a dielectric surface. For molecules, in contrast to atoms, the energy 
levels are closely separated and may thus easily allow transitions. For recent 
investigations along these lines, cf. Refs. [HI US].) With this, Eq. (JTJ) will 
be modified to Aa n = 'Yjm^ n a mb nm - Now, the state n does not only receive 
contributions, but gives away contributions to other states also. The latter 
must follow from the corresponding increase of probabilities for the other 
states. Omitting the latter for the moment, the perturbed coefficients are 

a ln = a n + Aa n = a n + ^ a m b nm . (11) 

m^n 

The a n will have complex phase factors, and in thermal equilibrium one must 
assume the phases of a n and a m (m ^ n) to be uncorrected. Thus by thermal 
average 

Ka m ) = 0. (12) 
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With this the new probability of the state n becomes 



Pin (^ln^l") | | ^ \@"m\ B nm , 

m^n 

B n m brim b nm I brim | • ( 1 3) 

The last term is the increase in probability from the other states. Likewise, 
the state n must obey a similar loss of probability to other states to conserve 
probability The loss to other states is thus J2 m ^n \ a n\ 2 B mn . With Eq. (jSJ) 
we have b mn = 6* m , by which B mn = B nm . The latter equation reflects that 
the transition probabilities between each pair of states are the same in either 
direction. With this, the resulting perturbed probability of state n modifies 
the expression (fT3|) into 

Pin = Wn\ 2 + $^(M 2 - Wn\ 2 )B nm = P n + J2( P m ~ Pn)B nm . (14) 

mj^n m 

The change in energy can now be evaluated as 

AE = E n(Pln ~Pn) = Yl E ^ R <n ~ Pn)B nm 
n nm 

^ ^ {.Bn Bm)PmB nm ~~\~ ^ ^ j yBmPm BnPn) B nm ^ Bm) PmB n m- 

nm nm nm 

(15) 

Utilizing the symmetry with respect to n and m in this expression and in- 
serting for P m from Eq. (jSj) we find 



AE = ^Yj{E n -E m ){P m -P n )B nm = ^J^^^^^nmSinhil^n^Br, 

(16) 



2 / J y ii lll/\ in ii J mil ^ 

nm nm 

with A nm = E n — E m , and where from Eqs. ffTUl) and ffT3"|) 

B nm = ^A nm An m q{—uj nm )q(u nm ). (17) 

Here it is to be noted that AE > 0. We conclude that whenever a system in 
thermal equilibrium is disturbed by some external perturbation the energy 
always increases (or, is unchanged), i.e., energy is dissipated. The dissipa- 
tion occurs to second order in the perturbation. To first order there is no 
dissipation; the changes are adiabatic. 
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We note in passing the similar nature of the formalism for an electromag- 
netic field in a dissipative medium: the mean quantity of heat developed per 
unit time and volume is 

Q = Lue"(u)(E 2 ), 

where e" denotes the imaginary part of the permittivity, e = e' + is" . Ir- 
reversibility of the dissipation process implies the condition e"(u>) > for 
positive cu. Cf., for instance, ref. [2D], Sec. 80. 



3 Energy dissipation from friction force 

In Ref. [13] , we evaluated the dissipated energy by considering the friction 
force. The q(t) denotes (or can be interpreted to denote) the position 

x = x{t) = q(t), (18) 

and A becomes correspondingly the operator for a force. By use of the 
Kubo relation [151 EE EH] f° r this situation the resulting force due to the 
perturbation is 

F f = I <t>A A {t-t')q{t')dt', (19) 



where 

0AA(t) = ^Tr{p[A,A(t)]}. (20) 



Here 



e 



-/3H 



p = with Z = Tr(e-P R ), 

Zj 



is the canonical density matrix, and 

A(t)=e itH/h Ae- itH/h , (21) 
with H the Hamiltonian. With velocity 

v(t) = x(t) = q(t) (22) 
the total energy dissipated by the system is 



/oo /* oc r pt 

v(t)F f dt = - / q(t)(f)AA(t - t')q(t')dt 

oo J -co J —oo 



dt, (23) 
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which is the same result as in Eq. (27) in Ref. [H]. Now the quantity q(t) 
need not be a position as given by Eq. (TTgj) . but as it can be interpreted as a 
position in a generalized sense, we concluded in [13] that the result (123)) has 
a broader applicability We will now show that this is actually the case, by 
showing that the result (123]) is the same as the new result (fT6l) obtained in 
the present work, by means of time-dependent perturbation theory. 
With wave function representation we first have 

e-? H ^J2Mx)e^ En r n (xi), (24) 

n 

pAA{t) = \Y,\ M*>~^Mxi)W m (xi)e™ mt fa(x 2 )A 

nmk 

x Mx2)e- iu,kt ^* k {x 3 )d Xl dx 2 . (25) 

Thus we obtain 

Tr(pAA(t)) = \Y. ^'"A.m^^e-^, (26) 

nm 

as J ipl(x)ip n (x)dx = 5kn (^3 — x\ — x), and A nm is given by Eq. (Q. 
Likewise we calculate Tr(pA(t)A) by exchange of u n and u m in Eq. ( |26l) . 
The response function becomes 

<t>AA{t) = ^Tr {p[A, A(t)}} = \Y1 M nm (e- lUnmt - e iawt ), (27) 



n in 



with 

M nm = -L e -¥iz n +E m ) sinh ^f3A nm )A nm A* nm (28) 

(recall that A nm = E n -E m = fiw nm , A mn = A* m ). The expression for M nm 
follows if one first exchanges n and m in Eq. (|26j) . then adds the resulting 
term to it and divides by 2. By inserting Eq. (12T|) into Eq. ( 12 3 p one gets the 
integral 

J = jf" y g(t)g(t') ^-^'e^' - e^e - *^) rft'c/t, (29) 

where here u = u nm [(E n — E m )/h = A nm /h}. By partial integration and 
exchange of integration variables t and t' (in the last term below ) we get 

/ = ico J J q{t)q{t') (e^V"*' + e iuJt e- luJt '^j dt'dt 



= IU 



/ / q(t)q(t')e- iujt e iuit 'dt'dt = iuq(u)q(-u). (30) 



J — oo J — oo 



By inserting this into Eq. (j2"3"j) via Eq. ( 1271) we get for the dissipated energy 



With u = co nm = A nm /h and M nm given by the expression (1281) this is 
nothing but the result (1TB]) together with (1171) obtained by time- dependent 
perturbation theory. Thus we have been able to derive the same expression 
for the dissipated energy in two independent ways. 

4 Friction between harmonic oscillators 

In Ref . [H] the friction between a pair of harmonic oscillators with interaction 
as in Eq. (J2J) (t > 0) was evaluated. Here, we will evaluate the energy 
dissipation by a direct use of Eq. (1161) or (13T1) . The first term in Eq. (T5]) 
can be disregarded as it gives a reversible force, distinct from dissipation. 
Further, we replace t with te~ vt (rj — > 0) to make the interaction vanish as 
t — > oo. 

For harmonic oscillators one can introduce the usual annihilation and 
creation operators 




(31) 



run 




(32) 



{i = 1, 2), with properties 



a 



n) = \/n + 1 \n+ 1), 



a\n) = y/n \n — 1). 



(33) 



With this the interaction becomes 



Aq(t) = 7(aict2 + dia 2 + a\a2 + alal) t e 



(34) 



where 



7 = (^Dh) 



1/2 (vW), ^ 



(35) 



2mim2o;ia;2 
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Since here only small r\ (— > 0) is considered, the terms aia 2 and a\a\ will not 
contribute. Thus we can use 

A = dial + a\d2, and q(t) = 'yte" 11 . (36) 

For the matrix elements we then get 

Au.na.ni+l.na-l = (nin 2 |oi4|ni + l,n 2 - 1) = a/«1 + 1 

-4m,n 2 ,ni-i,n 2 +i = (nin 2 |a{a 2 |ni - l,n 2 + 1) = y^V^ + 1, (37) 
while all other elements are zero. The Fourier transform of q(t) is, for t > 0, 

q(u) = 7 f° fe-*e-*"* = 7 (38) 

Jo iv + ™y 

so that, for 77—7-0, 

here w = U\ — u 2 , where U\ and u 2 are the eigenfrequencies of the two 
oscillators. Further, with uj — > (m = n ± 1) 

A nm sinh(l/3A nm ) -> = l -^±hw) 2 = (40) 



Then the matrix elements (1371) should be squared and averaged by the Boltz- 
mann distribution given by Eq. (fUJl) . We have (ni) ~ (n 2 ) ~ (n), with 
uj\ — > 0J2, and 



00 

cc \ ^ ~ X 



Z ^ 1-2 

n=0 



{in - — - > //./•" . x = e P hbJ1 



7 , 



z = ^Y,- n = Y^- x ( 41 ) 

n=0 

Then (n) + 1 = 1/(1 — x), by which 

((m + l)n 2 + ni(n 2 + 1)) = 2((n) + l)(n) 
2x 1 



;i-x) 2 2sinh 2 (±/^ i ; 
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(42) 



Inserting in Eqs. (IT5|) and ( 1P7|) by multiplying together Eqs. (I3"9"j) . (HO]) , and 
(I42p we obtain for the energy dissipation 



AE 



5(u)i - uj 2 ). 



(43) 



877 sinh 2 (|/3wi) 



With 7 inserted from Eq. (|35|) this is the same as the result (21) of Ref. [H] 
with its Eq. (19) for the friction force inserted. 

According to Eq. (j4"3~]) . the case of zero temperature (j3 — > 00) yields 
AS — > 0. This result is related to our assumptions, including slowly varying 
coupling or low velocities, i.e. 77 — > in Eqs. (155]) and (13"9"]1 . At more rapidly 
varying coupling or higher velocities also finite frequencies would contribute, 
leading to a finite energy change and a finite friction force at T = 0. 

5 Summary 

Let us summarize our work as follows (cf. also the end of Sect. 1): 

1. The total energy dissipation was calculated for a system perturbed 
by a time-dependent interaction. The change in energy is basically a second 
order effect but it was calculated with the use of standard time-dependent 
perturbation theory only, the reason being the absence of cross-terms due to 
uncorrelated phases of eigenstates. The energy change was always found to 
be positive or zero. The result agrees with, and confirms, our previous result 
of Ref. [H] obtained in a different and independent way. 

2. Our results are basically worked out at any temperature, and we 
have assumed initial thermal equilibrium. Moreover, we have assumed low 
velocities and nonrelativistic mechanics. Photons are accordingly not present 
in the theory. Photons were included, however, in our earlier study |16j . 

3. The energy change AE is finite in general. This corresponds to a 
finite friction force. In the limit T — > the expression ( 143]) gives, however, 
AE — > 0. This result is due to our assumption about constant velocity, 
involving slowly varying coupling. For couplings varying more rapidly, there 
will also be a friction force at T = 0, due to transitions to excited states. 
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